QED-based Optical Bloch Equations without electric dipole approximation: A model for a two-level 
atom interacting with a monochromatic X-ray laser beam 
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We derive a set of optical Bloch equations directly from the minimal-coupling Hamiltonian density of the 
bound-state quantum electrodynamics (bound-state QED). Such optical Bloch equations are beyond the former 
widely-used ones due to that there is no electric dipole approximation on the minimal-coupling Hamiltonian 
density of the bound-state QED. Then our optical Bloch equations can describe the time evolution of a two-level 
atom in a monochromatic light of arbitrary wavelength, which are suitable to study the spectroscopy and the 
Rabi oscillations of two-level atoms in X-ray laser beams since that the wavelength of X-ray is close to the atom 
to make the electric dipole approximation invalid. 
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I. INTRODUCTION 

Optical Bloch equations (OBEs) are the most widely used 
time-dependent equations for light-atom interaction in quan- 
tum optics til|2|], with which the stimulated absorption, stim- 
ulated emission, and spontaneous emission of light by atoms 
are well described, and the spectroscopy and Rabi problem of 
atoms can be studied well. The OBEs play the central role in 
both the semi-classical models and the quantum models such 
as Jaynes-Cummings model Q3D - In these models, the interac- 
tion Hamiltonian is H; = -<?D ■ E, which is restricted by the 
electric dipole approximation (EDA). Such approximation is 
a combination of two individual approximations. The first one 
is the dipole approximation, which requires that the radius of 
an atom r (~ O.lnm) is much less than the wavelength of the 
light A, 



2tt • - = k ■ r « 1 . 

A 



(1) 



The second approximation neglects the coupling between the 
atoms and the magnetic field of the light B due to that it is 
much smaller than the coupling between the atoms and the 
electric field E of the light. Therefore, these widely-used 
OBEs in quantum optics textbooks are only valid for the in- 
teraction between atoms and long-wavelength light. 

Recent realization of the free-electron lasers that work on 
X-ray wavelengths JH Ht] (although only femtosecond pulses 
are available at present) provides an opportunity for the exper- 
imental study the time evolution of atoms in monochromatic 
X-ray beam. With the steady-state solutions (f » r, with r 
being the lift-time of the excited state) and instantaneous so- 
lutions (t < t) of OBEs one can study the spectroscopy of 
atoms in continues-wave (CW) laser and the Rabi problem of 
atoms in ultrafast laser pulse (such as femtosecond pulses), re- 
spectively. Such OBEs are different from the X-ray scattering 
theories [6] due to that the OBEs are time-dependent which 
study the evolution of the system, while the X-ray scattering 
theories study the time-independent initial and final states of 
the system. However, the wavelength of X-ray (~ Q.Qlnm 
to 10«m) does not satisfy the dipole approximation condi- 
tion (r <K A), which makes the widely-used OBEs in quantum 



optics textbooks invalid in studying the interaction between 
atoms and X-ray laser. Then a new set of OBEs beyond the 
electric dipole approximation is required for the spectroscopy 
and Rabi oscillations of atoms in X-ray laser beams. 

Quantum electrodynamics (QED), which is the relativis- 
tic quantum field theory for the electromagnetic interaction, 
studies the light-matter interactions at the fundamental level 
B7|]. Bound-state QED is the extension of QED which studies 
the interaction between photons and bound electrons in atoms 
with the language of QED. The early bound-state QED mod- 
els are developed by Furry JH, Salpeter and Bethe JH Holl . 
as well as Gell-Mann and Low ll 111 to study the bound-state 
problem in quantum field theory. In the past years, many 
bound-state QED models were developed to correct the elec- 
trons' energy levels in atoms and ions IIL?T4l4tl . including few- 
electron atoms 11151 1 1 €h1 and many-electron atoms Ill7l - l22ll . 
Since the interaction Hamiltonian in the bound-state QED 
models does not have the EDA, it can describe the interaction 
between the bound electrons in atoms and the light at all wave- 
length. Then it is interesting to ask a question that whether a 
new set of OBEs can be derived directly form the interaction 
Hamiltonian in bound-state QED? If so, such OBEs can get 
rid of the EDA in order to describe the interaction between 
atoms and light of arbitrary wavelength, especially to meet 
the recent development of X-ray laser. 

In this paper, we start from the minimal-coupling Hamil- 
tonian density between the bound-state Dirac field and the 
free electromagnetic field in bound-state QED to derive the 
new style of OBEs. Our new OBEs can describe a two- 
level atom interacting with a monochromatic light of arbitrary 
wavelength, which meet the requirement of studying the spec- 
troscopy and the Rabi oscillations of two-level atoms in X-ray 
laser beams. Our paper is arranged as follows: In Section 
II, we present the interaction Hamiltonian for a light-atom in- 
teracting system form the interaction Hamiltonian density be- 
tween bound-state Dirac field and free electromagnetic field 
in bound-state QED, and derive the time-dependent equations 
of the ground and excited states of a two-level atom in inter- 
action picture. In Section III, we involve the density matrix 
of the two-level atom and the spontaneous emission of its ex- 
cited state into the time-dependent equations in last section, 



2 



and derive the OBEs for a two-level atom interacting with 
monochromatic light. We present how the new OBEs can 
describe the spectroscopy and the Rabi oscillations of a two- 
level atom in X-ray laser. Section V is the conclusion of this 
paper, which also gives the discussions on possible extensions 
of such new OBEs. 



II. TIME-DEPENDENT EQUATIONS FOR A TWO-LEVEL 
ATOM IN A MONOCHROMATIC LIGHT FORM 
BOUND-STATE QED 

In the interaction picture, the interaction Hamiltonian 
density of quantum electrodynamics (QED), which is the 
minimal-coupling between free Dirac field and free electro- 
magnetic field, is 



tii = -Li 



-eipYA^ip, 



(2) 



with <p being the Dirac field, y being the Dirac matrices, 
being the covariant four-potential of the electromagnetic field, 
and e being the elementary charge. The minimal-coupling 
Hamiltonian density for the interaction between the bound- 
state Dirac field and free electromagnetic field in bound-state 
QED has the similar form H[H 



tif = -eipeA^B, 



(3) 



Here <ps is the bound-state Dirac field which describes the 
bound electrons in atoms, and A 1 * is the four-potential of the 
free electromagnetic field. 

Quantization of the bound-state Dirac field, <fs, is different 
to it of the free Dirac field, tp. The reason is that the eigenstates 
of the bound electrons are the discrete wave functions of en- 
ergy and angular momentum rather than the continues wave- 
functions of kinetic energy and momentum for free electrons . 
Since no positrons exist in atoms, we simplify the Furry pic- 
ture flSQ by removing the creation and annihilation operators 
of positrons, and only preserving the creation and annihila- 
tion operators of electrons. Then the bound-state Dirac field 
for the valence electron in a two-level atom can be canonically 
quantized as 



<fB = Z Z Pp,m<l>m(x)e* px e* Emt b m , 



p m=e,g 



= Z Z P p „<r m {x)e-lr'e-i E ->bl 



(4) 



p m=e,£ 



Here, <p,„{x) is the space-dependent wave-function of the 
bound electron in the state \m > ( with \e >= \\ e , g > be- 
ing the excited state and \g >- \Q e , l g > being the ground 
state ), p is the external linear momentum of the bound-state 
Dirac field, E n = Hoj m is the total energy of the bound-state 
Dirac field in the state m, and p p<m is the normalization param- 
eter, bjjb,,, is the creation/annihilation operator of the bound 
electrons in the state \m >, which obeys the anti-commutation 
relation {b] n ,b m ,} = 8 mm! , and \bl„b],} = [b m ,b m >] = 0. 



The canonical quantized monochromatic electromagnetic 
field in Coulomb gauge can be written 



,1=1 



t -ik-x-o>,t 



+ ae' 



ik-x+tait 



), 



(5) 



with k being its wave-vector, <±>\ being its frequency, s = 
(2?^) 1 being the normalization coefficient, being the 
vector polarization, and eft la being the creation/annihilation 
operator of photons which obeys the commutation relation 
[a, a'] = 1, and [eft, a J ] = [a, a] = 0. Since the angular mo- 
mentum of a bound electron is fixed at a certain energy level in 
atoms, we assume that the polarization of the electromagnetic 
field (circular or linear polarization) just fits the requirement 
of the angular momentum difference {cr ± or n transitions) be- 
tween \e > and \g >. Then Xu=i §M can be replaced by eu for 
such two-level system. 

After the quantization of both bound-state Dirac field and 
monochromatic electromagnetic field, we can integrate Eq. Q 
in the three-dimensional space to obtain the time-dependent 
interaction Hamiltonian ///(f)- 

Hi(t) = j dx 3 ip B A fl ip B 

= -espipgablbge-^'-"*-"* 

x J f e {xTe k ( Pg < < x)e^-^ x+mx dx i (6) 

-eE^ e a%\b e e^-^-^ 

X j ^{x)h<pe{x)e^-^ )x - ikx dx\ 

Here we adopt the rotating-wave approximation to ignore 
the counter-rotating frequency parts exp [+i(oj e - co g + a>i)t], 
which is an widely used rational approximation in quan- 
tum optics 01 0]. The integrations of J (/>* e (x)ei l (p g (x)dx 3 and 
J (f> g (x)ek<pe(x)dx 3 , which show the transition probability be- 
tween \e > and \g >, can be calculated with the wave-function, 
(f>e(x) and 4> g {x). Here we define 

A — a>i — ((±> e — u) g ), 

q = -es/rj3 g j tiWerfgWefc'-^ ^dx 3 , (?) 

<t = -es/3* g /3 e J f g {x)h4>e{x)e^-P-> x - ik x dx\ 

The the spatial integrations (q and q*) are the key differences 
between the conventional interaction Hamiltonian with EDA 
and the minimal coupling Hamiltonian in Eq. (O. In the inter- 
action Hamiltonians which are based on EDA (HQ], the spa- 
tial integrations (q and q*) are not considered due to that the 
wavefunctions of electrons are localized to the spatial point 
where the atom is (x = 0), then such integration are re- 
placed by localized electric-dipole transition matrix elements 
(< e\r\g > and < g|r|e >) at x = 0. In the minimal coupling 
Hamiltonian, Eq. (|6), the integrations can be directly calcu- 
lated with the wavefunctions of electrons, which is free from 
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FIG. 1: (Color online) Transitions of a two-level atom in a quan- 
tized monochromatic light. Solid lines with arrows denote the stim- 
ulated excitation and emission of a monochromatic photon between 
the ground state |0 e , l s ,n+l > and the excited state |l e , g ,n >, where 
e and g are the excited and ground state, respectively, n is the parti- 
cle number state of the photons, u>i is the frequency of the photons, 
and A = a>i - (a) e - a> g ) is the detuning of the photons to the energy 
difference between |0 e , 1^ > and |l c , g >. Dotted line with an arrow 
denote the spontaneous emission of the excited state, with T being 
the spontaneous emission rate 



The set of two equations in Eq. ( fTTT i is an ideal model where 
the spontaneous emission of the excited state is not consid- 
ered. For a real system of a two-level atom interacting with 
monochromatic light, the spontaneous emission of the excited 
state cannot be ignored. The density matrix of the system is 
necessary to combine Eq. ( fTTT i and the spontaneous emission 
into a set of equations, which are just the optical Bloch equa- 
tions (OBEs). 



III. OPTICAL BLOCH EQUATIONS FOR A TWO-LEVEL 
SYSTEM BEYOND THE ELECTRIC DIPOLE 
APPROXIMATION 



To derive the OBEs, we need to introduce the density matrix 
for the two-level atom interacting with monochromatic light at 
first. The density matrix is 



the wavelength of light. Then Eq. (|6]l can be used to derive the 
OBEs with arbitrary light wavelength. With Eq. (|7) the time- 
dependent interaction Hamiltonian ///(f) can be rewritten as 



///(f) = q*e- ,A 'a ] 'blb e + qe ,A 'ab' e b g . 



(8) 



This is the bound-state QED based interaction Hamiltonian 
for a two-level atom interacting with photons. 

Following the standard procedure in quantum optics text- 
books [1, 2], we can define the time-dependent particle num- 
ber state |*P(f) > of the system in the interaction picture, 

\¥(t) >= c e (t)\l e ,0 g ,n > +c g (t)\0 e , l g ,(n + 1) >, (9) 

where c e (t) is the probability amplitude of the atom in the ex- 
cited state |l c ,0 ? >, c g (t) is the probability amplitude of the 
atom in the ground state |0 e , \ g >, and n is the number of the 
photons. The Shrodinger-like equations for ///(f) and \*¥(f) > 
in the interaction picture is 



ih-mt) >= Hjitm) > 

at 



(10) 



Then with Eq. (8), ®, and ( [Tol l, we can obtain the time- 
dependent equations of c e (t) and c g (t), 



c e (t) = --q Vn + le iA 'c g (t), 
c g (t) = --q* Vn + le~ iA 'c e (t)- 



(11) 



This coupled set of equations describes the Rabi oscillation 
between the two level, \l e ,0 g ,n > and |0 C , l g ,(n + 1) >, with 
the Rabi frequency 



Q 



2<? Vn + 1 



, 2q* + 1 
Q = — . 



(12) 



P = 



Pee Peg \ _ 
Pgr P 



c*(t)c e (t) c* e (t)c g (t) 
c*Jt)c e (t) c*(t)c g (t) 



(13) 



with p ee = c*(f)c e (f) being the probability of the atom in the 
excited state \e >= \l e ,0 g >, p gg = c* g (t)c g (t) being the prob- 
ability of the atom in the ground state \g >- \Q e , l g >, and 
Peg - P* ge = c t(t)c g (t) being the off-diagonal elements which 
determine the transition between \g > and \e >. An obvious 
relation is that p ee + p gg = 1 . 

After introducing the density matrix, we need to include the 
spontaneous emission of the excited state \e >- \ l e , Q g >. The 
standard theory of spontaneous emission is the Weisskopf- 
Wigner theory Il23ll . which presents that the spontaneous emis- 
sion is the result of the interaction between the two level 
atom at the excited state \e > and the vacuum state of multi- 
mode electromagnetic field. The former interaction Hamilto- 
nian in Weisskopf-Wigner theory is also based on the electric 
dipole approximation (EDA). In principle, we can rewrite the 
Weisskopf-Wigner theory with the minimal-coupling Hamil- 
tonian of the bound-state QED 



H,(t) = ^(qle 



■'a'b'gb, 



+ q k e lAt ab\b i 



(14) 



and the initial state of the system 

mi) >= c e (t)\l e ,0 g ,0 > + 2^(010., l g , h > ■ (15) 



Such new form of Weisskopf-Wigner theory can get rid of 
the EDA. However, this work is worth to do in the future and 
we do not to present it here. In this paper, we adopt a con- 
stant T as the spontaneous emission rate of the excited state 
\e >= \l e ,0 g ,n > (then the lifetime of the excited state be- 
comes 1 /F), which is the widely-used approach in quantum 
optics ll|,|2t]. 

Now with Eq. ( fTTT i and F, we can write the time-evolution 



4 



equations for every element of the density matrix, which are 



pee(t) = C* e (i)C e (t) + C* e (f)C e (f) - T Pe 



-q V« + lc 



Pge 



-a Vn + 1 
n 



e' 'p eg - Tp e , 



Pgg(t) = c*(t)c g (t) + c*{t)c g (f) + Tp e . 



— q Vn + 1<? 



iAt 



Peg + -Zi 



Vn + r< 



-iAt , 



+ rp 6 



p e ,(0 = (c* e (t)c g (t) + c* e (t)c g (t))e il 



+ (zA - -)p es 



= 7? 



lp^(0]* 



P 

Vn + l (Pgg ~ Pee) + OA - - 



)P> 



= l -q* Vn + l(p e , 



■ Pgg) ' 



eg ' 



(zA + -)p ge - 



(16) 

Here p e? = p eg e' A ', and p ?( , = p ge e ,At . The equations in 
Eq. ( [ToT l are our bound-state QED based OBEs for a two- 
level atom interacting with monochromatic light of arbitrary 
wavelength. We obtain these equations from the interaction 
Hamiltonian, Eq. ([8]), which is directly derived from the in- 
teraction Hamiltonian density between a bound electron and 
a monochromatic electromagnetic field in bound-state QED. 
Therefore, Eq. (fTST i is the set of OBEs beyond the EDA. 
Comparing with the former OBEs that derived from Jaynes- 
Cummings model yfl, the brief change of our OBEs is that 
the Rabi frequency under the EDA is replaced by the iq/h, 
where q is the spatial integration of all space-dependent wave- 
functions of both photons and bound electrons (see Eq. (0). 
This is the reason why our bound-state QED based OBEs are 
free from any spatial approximation, including the EDA. 

A very useful solution of Eq.tfToTi is the steady-state solution 
of p ee under the condition p ee (t) = p gg (t) = p eg (t) = p ge (t) = 0. 
The steady-state solutions of p ee and p gg are 



Pee = 



A\q\\n + l)/(h 2 F 2 ) 



1 + 8|<?| 2 (n + l)/(zi 2 r 2 ) + 4A 2 /r 2 ' 

Pgg — 1 — Pee 

_ 1 + 4\q\ 2 (n + l)/(/z 2 r 2 ) + 4A 2 /r 2 
~ 1 +8\q\ 2 (n+ l)/(/i 2 r 2 )+4A 2 /r 2 ' 



(17) 



Fig. |2]shows the varying of p ee with A/T under two different 
photon numbers n and n' . The distribution of p ee is a Lorentz 
line-shape, which determines the amplitude and width of the 
absorption signal of the monochromatic light by a two-level 
atom in laser spectroscopy. The resonance point of the light 
and the two-level atom is at A = 0, where p ee has its maximum 



value 



4M 2 (»+i) 



rrfo ol .,. — r-. We can see such maximum value of p ee is 
always less than 0.5 due to that the n — > oo limit of p ee is just 
0.5 at A = . This means that the population of the excited 
state p e e is always smaller than it of the ground state p gg in 
a steady-state two-level atom. Besides, n and |g| 2 also deter- 
mine the full-width-at-half-maximum (FWHM) of the absorp- 
tion signal. In the condition \q\ 2 (n + 1) <K h 2 T 2 , the FWHM 
of the absorption signal is approximately T, which means that 
the absorption signal of a weak-intensity light has the natu- 
ral line-width F of the excited state. When the light intensity 
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FIG. 2: (Color online) The variation of the excited state p ee of the 
two-level atom with A/T in the steady-state limit (t » 1/F), with A 
being the detuning of the light and T being the spontaneous emission 
rate. The solid (black) line is for photon number n = Ti 2 T 2 /4\q\ 2 — 1, 
and the dashed (red) line is for the photon number n' = 2h 2 r 2 /\q\ 2 - 
1. Since the intensity of the monochromatic light is proportional to 
its photon number n in the particle number representation, a larger 
n leads to a larger maximum value and a larger full-width-at-half- 
maximum (FWHM) of the p ee 



(proportional to n + 1) increases, the FWHM of the absorption 
signal is broaden to the value -\/4\q\ 2 (n + l)//z 2 + F 2 , which is 
called power broadening in laser spectroscopy. Such steady- 
state solution of pee is also proportional to the radiation force 
F = fiTpee from the monochromatic light to the atom, which 
is very important in the laser cooling of atoms. 

Another useful solution of Eq. ( [ToT l is the evolution of p ee (t) 
before the system reaches its steady-state. From Fig.|3]we can 
see such evolution of p ee (t) is also determined by the values 
of n and \q\ 2 . When \q\ 2 (n + 1) = ft 2 T 2 /4, the p ee {i) reaches 
its steady value p ee ~ 3.3 near t = bT, and has less than one 
oscillation period before it; when \q\ 2 {n + 1) = H 2 F 2 , the p ee (t) 
reaches its steady value p ee ~ 4.4 near t - IT, and has two 
Rabi oscillation periods before it; when |g| 2 (n+l) = 47z 2 F 2 , the 
Pee(t) reaches its steady value p ee ~ 4.8 near t = 8F, and has 
five obvious oscillation periods before it; when \q\ 2 {n + 1) = 
16S 2 F 2 , the pee(t) reaches its steady value p ee ~ 5.0 near 
t = 9r, and has ten obvious Rabi oscillation periods before 
it. The larger values of \q\ 2 (n + 1) cause (i) larger steady-state 
values of p ee (t — > oo), (ii) larger damping time of p ee (t) before 
reaching steady-state, (iii) and larger Rabi oscillation frequen- 
cies as well as amplitudes of p ee (t) during the damping time. 
The value of \q\ 2 is determined by the space-dependent wave- 
functions of the excited and ground states (see Eq. (Q). When 
the Rabi oscillation disappears, the system reaches its final 
steady-state, where the phase parameter, e ±lAt , disappears and 
the dephasing effect due to the spontaneous emission happens 
ll24ll . Therefore, the way to increase the Rabi oscillation fre- 
quencies and amplitudes, as well as the damping time of p ee (t) 
is to increase the photon number «, which is just to increase 
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FIG. 3: (Color online) Evolution of the excited state population of 
p ee (t) between t = and / = 9/T under different values of (n + l)\q\ 2 
(with n being the number photons which is proportional to the inten- 
sity of the light). 1/r is the lifetime of the excited state |l e ,0 g , n >. 
A = 0. The dotted (black) line, dashed (red) line, dash dotted 
(green) line, and solid (blue) line are the evolution of p ee (f) with 
\q\ 2 (n + 1) = n 2 T 2 /4, \q\ 2 (n + 1) = TrT 2 , \q\ 2 (n + 1) = Ah 2 T 2 , and 
|<?| 2 (n+ 1) = 16/rT 2 , respectively. A larger |g| 2 («+ 1) leads to a larger 
steady-state value of p ee (t — > oo), larger Rabi oscillation frequency, 
and larger amplitude of p ee (i) before t = 9 IT. 



the intensity of the monochromatic light in practise. 

Our new style of Optical Bloch Equations reach the simi- 
lar conclusions to that of the Jaynes-Cummings model |3|] in 
the case when the wavelength of the light is much larger than 
the radius of an atom (ultraviolet light, visible light, infrared 
light, and microwaves). In other cases when the wavelength of 
the light is close to or smaller than the radius of an atom (X- 
rays and Gamma rays), the Jaynes-Cummings model is not 
suitable due to that the EDA becomes invalid. However, our 
new OBEs can still reach the results of Fig.|2]and Fig.|3]with- 
out the EDA. Then our OBEs can be applied to describe the 
evolution of a two-level atom (with transition frequency be- 
ing in the X-ray range) in a monochromatic in a X-ray laser 
beam. Such two-level atom could have the absorption spec- 
troscopy of X-ray laser that described by Fig. |2) in the steady- 
state (f » 1/E), and have the Rabi oscillation that shown in 
Fig.|3]before the two-level system reaching its steady-state. 

Another useful result in Fig. [3] is the population inversion 
(p ee > 0.5) of the two-level system during the first period of 
Rabi oscillations (t «: 1 /F). Since the Rabi frequency and am- 
plitude are proportional to V« + 1, the time from t — to the 
maximum value of population inversion (we define it as At) 
can be shorten by increasing the intensity of the monochro- 
matic light, while the p ee is also increased. Then the value of 
At can be obtained from the evolution of p ee {t). This value 
is useful for the pulse pumped atomic X-ray laser system 
JH]. For example, a three-level system with two ground states 
\g\ >> \gi >> and one excited state \e > is prepared with all the 
papulation being at \gi >. We can pump the population from 
the |gi > to \e > to make a population inversion between \e > 



and \g2 >, then the X-ray laser with the frequency between 
[e > and \g2 > can be made. In the case when the spontaneous 
emission rate from \e > to \gi >, Fi, is much larger than it 
from \e > to \g2 >, T2, the transition between \e > to \g\ > by 
a free-electron laser can be considered as a two-level system 
approximately, and can be described by the OBEs in Eq. ([ToT i. 
When the duration of an intense pulse from the free-electron 
laser is close to At, the population inversion between \e > to 
|gi > happens (green line in Fig. |3J, which is the most effec- 
tive population inversion between \e > and \g2 >■ Then from 
the OBEs in Eq. ( fl6] i one can calculate out At to optimize the 
duration of the pump pulse in order to give an efficient pump 
rate for the atomic X-ray laser. 

IV. DISCUSSION 

We consider a simple case that a two-level atom interact- 
ing with monochromatic light in this paper. Our work can be 
extend to the case of a multi-level atom interacting with multi- 
frequency light, with which the non-linear effect of the light- 
atom interacting system can be studied. The quantization of 
multi-frequency electromagnetic field is 

2 

A " = Z Z + ( 1 8) 

k A=l 

where tok = ck. The quantization of the bound-state Dirac 
fields in a multi-level atom is 

M 

ifiB = 2] Y J /3p,m4 > m(x)e^ px eT- E -'b m , 

(19) 

p m= 1 

where m is the subscript of a energy level, and M is the num- 
ber of energy levels. An example is the model of a three-level 
atom (\g >,\el >,\e2 >) interacting with a two-frequency 
(<D[,(l>2) light, with co\ (102) light causing the transition be- 
tween \g > and \el > (\e2 >). In such case, the time-dependent 
interaction Hamiltonian in Eq. ^ can be rewritten as 

Hj(t) = q y e iklt aib\j3 g + q 2 e iA2 'a 2 bl 2 b g + H.c, (20) 

where Ai = a>\ — {cj e \ — io g ), A2 — a>2 — (co e 2 — w g ), and q\ 
(^2) is the spatial integration in Eq. with /3* and 4>* e (x) be- 
ing replaced by j3" el (J}* e2 ) and f el (x) (fjx)). After defining 
the three-level density matrix and the spontaneous emission 
rate of |el > and |e2 >, which are Ti and F2 respectively, a set 
of OBEs can be derived from the interaction Hamiltonian in 
Eq.[20]to describe a typical nonlinear spectrum effect, which is 
the electromagnetically-induced transparency (EIT) ll25ll . Be- 
sides, the extensions of our QED-based OBEs with different 
number of atomic energy levels and light frequency can study 
other nonlinear spectra of atoms including recoil-induced res- 
onances (RI R) 12611 and electromagnetically-induced absorp- 
tion (EIA) [27]. Since all these nonlinear spectra are for- 
merly studied by the OBEs with EDA, when EDA becomes 
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invalid the extension of our QED-based OBEs would be a bet- 
ter choice. 

The QED-based OBEs also fully describe the interactions 
between bound electrons and both the electric and magnetic 
field of light. Besides the dipole approximation, the electric 
dipole approximation (EDA) also ignores the interaction be- 
tween the atoms and the magnetic field B of the light due to 
that this interaction is much smaller than it between the atoms 
and the electric field E of the light. However, the interac- 
tion Hamiltonian density in Eq. ([3]) is the minimal-coupling in 
bound-state QED, where the bound-state Dirac field <p B cou- 
ples the four-potential of the free electromagnetic field. 
The four-potential A M = A + (p contents the contributions from 
both the electric field of the light by E = -V0 - (1/c) ■ dA/dt, 
and the magnetic field of the light by B = V x A. Therefore, 
the bound-state QED based OBEs are beyond the EDA not 
only by the dipole approximation part, but also by the electric 
approximation part. 

In summary, we derive a new set of optical Bloch equa- 
tions (OBE) from the minimal-coupling Hamiltonian density 
between the bound-state Dirac field and the free electromag- 
netic field in bound-state QED, which is Eq. (f3). The time- 
dependent interaction Hamiltonian, Eq. ©, is the spatial inte- 
gration of the minimal-coupling Hamiltonian density, which 
does not have the electric dipole approximation (EDA). Then 
our new set of OBEs is an useful tool to study the quantum op- 



tics of the light-atom interacting system where the wavelength 
of the light is close to or smaller than the scale of an atom. 
Since recent realization of X-ray lasers [41 5] provide the op- 
portunity to study the quantum optics phenomena of atoms 
interacting with X-ray experimentally, our new OBEs can pre- 
dict that a two-level atom with X-ray transition frequency in 
a X-ray beam will have the absorption spectra which are de- 
scribed by Fig. [2] and the Rabi oscillation which are described 
Fig. [3] The duration of the free-electron laser pump pulse for 
the atomic X-ray laser can also be optimized by the solutions 
of the new OBEs. When the new OBEs are extended to the 
case of a multi-level atom interacting with multi-frequency 
light, the non-linear effect of such light-atom interacting sys- 
tem can be studied by them. Therefore, we hope our work 
could be an opening work of developing more OBEs beyond 
the EDA. 
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